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We report on a numerical study intended to examine the possibility that magnetic oscillations persist in 
type II superconductors beyond the point where the pairing self-energy exceeds the normal state Landau 
level separation. Our work is based on the self-consistent numerical solution for model superconductors of 
the Bogoliubov-deGennes equations for the vortex lattice state. In the regime where the pairing self-energy 
is smaller than the cyclotron energy, magnetic oscillations resulting from Landau level quantization are 
suppressed by the broadening of quasiparticle Landau levels due to the non-uniform order parameter of the 
vortex lattice state, and by splittings of the quasiparticle bands. Plausible arguments that the latter effect 
can lead to a sign change of the fundamental harmonic of the magnetic oscillations when the pairing self- 
energy is comparable to the cyclotron energy are shown to be flawed. Our calculations indicate that magnetic 
oscillations are strongly suppressed once the pairing self-energy exceeds the Landau level separation. 
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I. INTRODUCTION 

de Haas- van Alphen (dHvA) oscillations in, the mixed 
state of type-II superconductors, discoveredEl in NbSe% 
some time ago, have recentlyBij been observed in several 
additional materials. The oscillations are damped rela- 
tive to those in the normal state and become unobserv- 
able at sufficiently weak external magnetic fields. These 
findings have led to a number of theoretical studies of 
the modification of-naniial state Landau level structure 
in the mixed state.ETt3 Conclusions from these studies 
are not always completely consistent and no widely ac- 
cepted picture which covers all regimes of magnetic field 
has emerged from this work. Recently we reported on 
a thorough numerical study of the quasiparticle band 
structure obtained by solving the Bogoliubov-de Gcnncs 
(BdG) mean-field equations in the vortex lattice state of 
a simple model two-dimensional superconductor.E3 We 
found that at fields near H C 2, magnetic oscillations were 
clearly present, but that these were rapidly damped as 
the superconducting self-energy strengthened at weaker 
fields. We argued and found partial numerical support 
for the assertion that the effect of superconductivity was 
similar to the effect of a disorder broadening of the nor- 
mal state Landau levels, proportional to the pairing self- 
energy times Tin ^ 4 where is the Landau level index 
at the Fermi level. We also found that once the pairing 
self-energy became comparable to the Landau level sep- 
aration, the quasiparticle electronic structure in the vor- 
tex lattice state entered a complicated crossover regime 
which simplified with increasing pairing self-energy only 
when unambiguous vortex cores with associated bound 
states emerged. While magnetic oscillations were essen- 
tially absent once the vortex cores became distinct, we 
were unable to draw any clear conclusions concerning 



magnetic oscillations in the crossover regime. These cal- 
culations did indicate the possibility of a phase shift of 
7r for magnetic oscillations in the crossover regime, but 
the origin of this phase shift was not understood. Such 
a phase shift was also found in earlier work by Maniv et 
so based on an expansion of the free energy to fourth 
order in the order parameter. Recently, Maniv et al.t3 
have attributed this phase shift to a splitting of Landau 
levels in the vortex lattice state which they associate with 
the two vortices per electron magnetic flux quanta in the 
vortex lattice state. This suggestion has motivated us to 
examine the crossover regime in greater detail. 

Our study is based on numerical solution of the BdG 
equations^ for a model superconductor with a BCS pair- 
ing interaction, i.e. a 5- function attractive interaction 
modified by an energy cut-off. We solve the BdG equa- 
tions in a Landau level basis so that the band energy 
quantization which is the source of magnetic oscillations 
is incorporated in an exact way. The formalism neces- 
sary to carry out these calculations-in-a convenient way 
is fully described in our earlier workE-3113 and briefly sum- 
marized below. This approach necessitates a number of 
practical limitations on the scope of our study, (i) Nu- 
merical problems which arise because of oscillations in 
high Landau index quantum wavcfunctions make it con- 
venient to restrict our attention to single-particle states 
with Landau level indices smaller than w 60. (ii) We 
consider only two-dimensional electron systems; adding a 
third dimension creates no formal difficulty but does add 
to an already considerable computational burden, (iii) 
We approximate the magnetic field by its spatial average. 
The most serious of these limitations is the restriction to 
moderately large Landau level indices. Two-dimensional 
models will, if anything, overestimate the importance of 
magnetic oscillations and are even appropriate for some 
systems of current interest. The screening corrections 
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to the uniform external magnetic field are small close to 
H C 2 and are approximately uniform themselves, except 
for external fields close to 

In the present study Zeeman splitting is ignored since 
its effects are well understood. Most of the results we 
discuss use the grand canonical ensemble rather than 
the canonical ensemble appropriate to experimental sys- 
tems, since this eliminates the problem of determining 
the chemical potential self-consistently. In the normal 
state, there is little difference between magnetic oscil- 
lations in canonical and grand canonical-ensembles for 
Landau level indices larger than about sixJ13 In the mixed 
state, however, canonical and grand canonical ensemble 
results may differ. We have therefore, in some cases, 
executed the Legendre transform from the grand canon- 
ical to the canonical ensemble numerically in order to 
quantify the importance of magnetic oscillations in the 
chemical potential. 

In Section II of this paper we summarize the BdG for- 
malism which is the basis of our numerical calculations. 
In Section III we discuss the devolution of the Landau 
level structure in the quasiparticle spectrum as the su- 
perconducting order strengthens. It is this devolution 
which underlies the damping of dHvA oscillations in the 
mixed state. We find that a picture in which the normal 
state Landau levels simply broaden captures little of the 
process, hence the substantial difficulty in developing a 
simple analytic theory for the influence of superconduc- 
tivity on dHvA oscillations analogous to the simple and 
successful theory for the influence of disorder. In this sec- 
tion we discuss a plausible approximation which suggests 
that magnetic oscillations in the vortex lattice state in 
the crossover field regime will differ by a sign from those 
in the normal state. The possible sign change is associ- 
ated with a splitting in the density of quasiparticle states 
associated with each Landau level at fields below H C 2- 
We explain the origin of this splitting and comment on 
the failure of the commonly used diagonal approximation 
for the quasiparticle spectrum. In Section IV we care- 
fully examine magnetic oscillations in this regime and 
find that the sign change does not survive a more thor- 
ough analysis. Instead, the fundamental harmonic of the 
magnetization is strongly damped. The magnetization 
in this regime has substantial variation with field but the 
indications from our numerical calculations is that the 
field dependence is aperiodic. We conclude in Section V 
with a brief summary. 



II. BOGOLIUBOV-DE GENNES FORMALISM 

In zero field BCS theory, the pairing self-energy cou- 
ples only single-particle states at wavevectors k and —k. 
The property that the coupled states have the same band 
energy is favorable for the formation of a condensate 
of electron pairs. In a magnetic field the loss of time- 
reversal invariance makes it impossible to achieve this 



situation. The center-of-mass momentum of a pair of 
electrons, which is zero for condensate pairs at zero mag- 
netic field, has quantum fluctuations in a magnetic field 
~ h/£ where I — (Tic/ eB) 1 / 2 is the quantum magnetic 
length. Associated quantum fluctuations in the momenta, 
of the individual electrons contributing to the pair leadE£l 
to pairing between electrons in different Landau levels 
and therefore with different single-particle energies. It 
is this qualitative difference which is responsible, from 
a microscopic point of view, for the decrease of T c in 
a magnetic field. The well known dependence of H C 2 
on field, obtained from semiclassical theory or (near T c q) 
from Ginzburg-Landau theory, reflects in the microscopic 
theory primarily contributions from pairing between elec- 
trons in different orbital Landau levels .Ej It is not pos- 
sible to understand the modification of Landau levels by 
superconductivity, even in the regime,-iiear H C 2, unless 
one includes these 'off-diagonal' terms.Ej 

The BdG mean-field equations for a superconductor 
in a constant magnetic fieldli3i23 replace the 2x2 sec- 
ular matrix of BCS theory at zero magnetic field by a 
secular matrix of order 2N (where N is the number of 
Landau levels within a pairing cut-off energy) for each 
wavevector, k, in the Brillouin zone of the vortex lat- 
tice. The diagonal (normal) electron and hole blocks of 
the secular matrix are diagonal in the Landau level basis 
with elements given by £„ and — respectively where 
£ n = (n + l/2)hu> c — fx. Here lo c = eB/mc is the cy- 
clotron frequency and /x is the chemical potential. It is 
the simplicity of the diagonal block which makes such a 
basis convenient—The off-diagonal (pairing) blocks have 
matrix elementstZl 
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In these equations A is the BCS coupling constant so that 
\fiLo c = V/(2irl 2 ) where V is the strength of the attrac- 
tive interaction. The vortex lattice primitive vectors are 
{0,a y ) and (a x , —a y /2) with a x a y — nl 2 (a x — \[%a y j2 
for a triangular lattice). The sum over j in Eq. 1 is over 
the possible partitionings of the total quantized kinetic 
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energy of the pair, hoj c (N + M + 1), into contributions 
from the pair center of mass motion, hu> c (j + 1/2), and 
the pair relative motion, Tllo c (N + M — j + 1/2), with 
{D^ M ) 2 the probability that a pair of electrons in Lan- 
dau levels N and M will have center-of-mass kinetic en- 
ergy 7kJ c (j + 1/2). 

In this formalism the order parameter in the vortex lat- 
tice state is paramaterized by a small set of numbers, Aj, 
which should be determined by solving the BdG equa- 
tions self consistently :EZl 

4 3 = -ErE|xw,-i(^) 

NM k 

where is the /x'th positive eigenvalue of the secular 
matrix, (u M -,v^ J) is the corresponding eigenvector, and 

' v Nk Mk' 1 ° ° ' 

/i 1 is the Fermi function. Nk = L x L y / (2irl 2 ) is the num- 
ber of k points {L~Ly is the area of the system). The 
Abrikosov solution^ for the order parameter near H C 2 
corresponds to a solution with only Ao 7^ and it is easy 
to verify that this solution is recovered in the appropri- 
ate limit. For a triangular flux lattice, the lowest energy 
solution has A 3 real and non-zero only for j — 6m where 
m is an integer. 

We determine the magnetization by numerically differ- 
entiating the appropriate thermodynamic potential with 
respect to magnetic field. The grand potential may be 
expressed in the following form which we use for our nu- 
merical calculations:t3 

fl = J2^NN N +E P -TS (6) 

N 

where the pairing self-energy is 

x 3 

Here Nm is the occupation number of Landau level N 

fj,k 

and S is the entropy: 

S = -^B 1 -/f)Ml-/f)+/£M. (9) 

fj,k 

For canonical ensemble calculations we calculate the free 
energy F over a range of electron densities from the grand 
potential calculated over a range of chemical potentials 
by using 

F = Q + vJ2 Nn - ( 10 ) 

N 



where both F and the density depend parametrically on 
/j,. The canonical ensemble magnetization is determined 
by numerically differentiating F with respect to field at 
fixed density. A portion of the discussion of our results 
is motivated by an equivalent alternate expression for 17 
in terms of quasiparticlc energies lEj'LrEj 

9k T E^ 

k fik B n 

The last term here is a double counting correction for the 
pair interaction energy. 

The magnetization is determined by numerically differ- 
entiation: M(B) — —dft/dB. In practice, we generate 
results as a function of = u/hcu c — 1/2 and calculate 
energies per state in the Landau level in units of Tiuj c . 
Therefore, the derivative for the magnetization has two 
terms, the first coming from differentiating an explicit 
dependence on B (that is, fl = ^IqB 2 , with one power of 
B coming from huj c and the other from the Landau level 
degeneracy factor), the second from the dependence of 
f2o on which is determined numerically. Note that 
we do not need to perform separate calculations to de- 
termine the density dependence of Q mentioned above 
and the field dependence of Q required for the magneti- 
zation. Similarly, in the canonical ensemble the magneti- 
zation can be expressed in terms of the derivative of the 
corresponding dimensionless free energy with respect to 
N = En n n- 

III. LANDAU LEVEL DEVOLUTION IN THE 
MIXED STATE 

We first analyze the secular matrix in the limit of small 
Ao. Our objective here is to understand the behavior 
of the mixed state quaisparticle bands over one period 
of the normal state magnetic oscillations. Consider the 
case where — n (n an integer). For this case, we 
note that for each electron energy in the upper diagonal 
block, there will be a hole energy of the same value in the 
lower diagonal block. For the Landau level at ji, these 
two have the same index (n), otherwise, their indices are 
different (n+m and n—m). When the order parameter is 
small, the strongest mixing of a particle in Landau level 
n + m will be with a hole in Landau level n — m. The 
degeneracy of the particle and hole levels will be lifted 
by the matrix elements in the pairing block which are, in 
general, off diagonal in Landau level index. At a given k 
the two levels will be split by 2|F n+mj „_ m | for all Lan- 
dau levels within the pairing cut-off. In particular, one of 
the quasiparticle energy levels at zero in the normal state 
will be shifted up by l-Fn.nl while one of the quasiparticle 
levels at huj c in the normal state will be shifted down 
by |F n +i )T ,_i|. Obviously, this splitting cannot continue 
to grow indefinitely since these two levels will eventually 
approach each other, leading to an avoided crossing. A 
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similar degeneracy occurs when = n + 1/2 with the 
electron level at Landau level index n + 1 + m and the 
hole level at Landau level index n — m being degenerate 
resulting in a similar splitting of each Landau level. In 
this case, one of the two quasiparticle levels which has 
energy l/27kJ c in the normal state will be shifted down 
by |F n +i, n |. For = n + 1/4 (or n + 3/4), the level re- 
pulsion effect is weak; that is, the Landau level splitting 
is most pronounced when degeneracies occur in the nor- 
mal state, partially invalidating the, analogy to Zeeman 
splitting suggested by Maniv et al$3 

We illustrate these points in an approximation where 
all matrix elements in the pairing blocks are taken to be 
the same constant, — l/2(n A1 7r) _1 / 4 Ao?iw c which is the 
large N limit of the matrix element At the chemical po- 
tential if A = 1 and x is se t to unity.Ea For cases consid- 
ered in this paper, we take the cut-off, wp, to be 1/2/i 
(thus for — 20, Landau levels 10-30 are involved in the 
pairing). The resulting eigenvalues for the above three 
cases are plotted as a function of Ao in Fig. 1. In the 
large A limit, oscillations in the low-energy quasiparti- 
cle eigenvalue spectrum are about the same magnitude 
and shifted by half a period relative to the oscillations in 
the normal state. The expression for the grand poten- 
tial in terms of quasiparticle energies suggests that this 
might lead to a n phase shift in the Fourier transform of 
the magnetization relative to the normal state case, i.e. 
to a change in sign of the oscillatory contribution to the 
grand potential. We say "might" since it is not obvious, 
even from Eq. 11, that a phase shift in the oscillations 
of low-energy quasiparticle energies will necessarily show 
up as a phase shift of the magnetization. (Eq. 11 in- 
volves three terms and each contributes strongly to the 
oscillatory dependence of M(B).) 

To examine this idea in more detail we have solved 
the BdG self-consistently at several different A values for 
rip G (20, 21) and calculated coefficients of the Fourier ex- 
pansions of quantities of interest within this interval. In 
the normal state the Fourier expansion coefficients vary 
slowly with the Landau level index associated with the 
interval over which the Fourier transform is performed, 
since the dominant variation with rip is periodic. In the 
mixed state we will have to check for this periodicity by 
verifying that the Fourier expansions in successive inter- 
vals are similar. We focus on the coefficient of the leading 
sine term in the Fourier expansion which is the dominant 
term in the normal state and refer to the Fourier expan- 
sion coefficients as harmonics of the magnetic oscillation; 
the terminology anticipates a repetition of the same pat- 
tern in successive intervals which does not always occur 
as we discuss in further detail below. For the interval 
rip G (20, 21) we find that the zero of the fundamen- 
tal sine harmonic of the Fourier transform of M(B) in 
this interval does indeed closely correspond to the point 
where the three curves in Fig. 1 cross. To test the degree 
of correspondence between the total oscillatory contribu- 
tion to the grand potential and the contribution from the 
lowest band of quasiparticles, we have also verified that 



the quantity 

k T El 
E x = — V ln[2 cosh( Ml (12) 

k 

(with 1 denoting the lowest quasiparticle band) has a 
magnetization whose fundamental sine harmonic agrees 
quite closely with that of the total magnetization. Note 
that Ei is essentially -1/2 the mean of the energies of the 
lowest quasiparticle band suggesting that there is some 
validity in associating magnetization oscillations with os- 
cillations in the low-lying quaisparticle bands. (This sim- 
ilarity of leading harmonics occurs even though the shape 
of the two 'magnetizations' with respect to rip are quite 
different; the correspondence does not hold for higher 
harmonics). Finally, we again note the qualitative dif- 
ference between Fig. 1 and what would be expected if 
the splittings were simply proportional to A as in^-jthe 
Zeeman-splitting analogy proposed by Maniv et al&3 In 
this case, avoided crossing effects at larger A do not oc- 
cur and additional zeroes would occur in the harmonics 
at larger A. 

The behavior seen in Fig. 1 should be contrasted with 
the commonly used diagonal approximation, where the 
only elements retained in the pairing blocks are diagonal 
in Landau level index. In this case, the eigenvalues are 
simply shifted from £/v to a/Cat + \Enn\ 2 - In this approx- 
imation the level splitting effect occurs only when £jy = 0; 
otherwise all quasiparticle Landau levels arc shifted away 
from the Fermi level. Because of this qualitative failure, 
we do not feel that the diagonal approximation is useful 
for understanding the electronic structure of the vortex 
lattice state except for the Landau level closest to the 
Fermi level and then only when n — rip. 

To examine how Fig. 1 is changed when the con- 
stant matrix element approximation is abandoned and 
details of pairing in the vortex lattice state are properly 
accounted for, we have solved Eqs. rL-|4 as a function of 
Ao. As discussed in our earlier work,t2lthe use of a sharp 
cut-off when solving the secular matrix leads to spurious 
effects in M(B) associated with-jkhe ratio of the cut-off 
energy to the cyclotron energy.E3 To eliminate this, we 
elect to use a smooth cut-off with the pairing interaction 
between Landau levels N and M scaled by VWmWm 
where 

W N = 1.55e' iiN/0 - 5uiD)i (13) 

In Fig. 2 we show a plot of the density of states for 
rip = 20 and AAo = 1. Each quasiparticle Landau level, 
not only the Landau level closest to the Fermi energy, is 
split into two roughly symmetric subbands. This split- 
ting is due to particle-hole mixing. We have been unable 
to uncover a detailed connection between this splitting 
and the fact, emphasized by Maniv et al., that two su- 
perconducting flux quanta pass through each area of the 
vortex lattice state enclosing one electronic flux quantum. 

In Fig. 3, we show results for the vortex lattice quasi- 
particle bands which are analogous to those of Fig. 1 
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obtained using the constant matrix element approxima- 
tion. The plotted eigenvalue in this case is the mean 
eigenvalue of the lowest band using a 66 k point grid 
in the irreducible triangle (1/12) of the vortex lattice 
magnetic Brillouin zone. The results look very similar 
to Fig. 1 up to the point where the curves cross. This 
crossing point is close to the point where the spatially av- 
eraged pairing self-energy in a coordinate representation 
(Fq ~ 0.44AA ?ia; c ) is equal to Huj c (F being the vor- 
tex lattice analogue of the BCS gap). For smaller values 
of AA , we are in the quantum regime where we expect 
strong magnetic oscillations. As AAo increases the de- 
pendence of the eigenvalues on n M weakens and magnetic 
oscillations are correspondingly damped. The oscillations 
are further damped in this regime by the non-zero width 
of the Landau levels which reflects the the non-uniformity 
of the order parameter. The width is linear in Ao for 
small Ao and should lead to an exponential suppression 
of magnetic .oscillations with an effective scattering rate 
linear in Ao Jl3 At higher values of AAo we initially enter 
into the crossover regime and then into the regime where 
well-defined vortex cores emerge. The fact that the mean 
eigenvalues increase in this regime reflects the crossover 
of the lowest ^energy quasiparticle states to vortex-core 
bound states.Ej The eigenvalues clearly still have a sub- 
stantial dependence on rt M within the interval (20, 21), at 
least in the crossover regime, although the dependence 
is much weaker than in the constant matrix element ap- 
proximation. 

Up to this point we have been performing calculations 
at fixed AAo, i.e. at fixed pairing self-energy. To com- 
pute the magnetization we should in principle determine 
Ao self-consistently at each value of n M and keep A fixed. 
To facilitate comparisons with the preceeding results for 
the quasiparticle bands we have chosen instead to allow 
A to vary with n M so that self-consistency is achieved 
at a fixed value of AAo. This self-consistent value of 
A at a fixed AA (A) is easily determinedE3 by using 
Eq.(5) to calculate the output value (Ag ut ) at A = 1: 
A = AA "/Ag ut . Results are shown in Fig. 4 for the 
fundamental sine harmonic of the Fourier transforms of 
M(B) and E\ versus AAo. A zero in the harmonic of 
M(B) occurs for Fq ~ 1.6Tiuj c (similar results are found 
for self-consistent calculations at fixed A). The zero of the 
harmonic of E\ is close to the zero for M(B) as claimed 
earlier. We note that in the small Ao regime, the depen- 
dence of the harmonic on Ao contains both linear and 
quadratic terms. The calculations are consistent with 
a crossover from quadratic to linear behavior when the 
quantity Fo/nj/ 4 exceeds 27t/cbT. We also present in Fig. 
4 canonical ensemble results generated from the grand 
canonical calculations by a Legendre transform. Devi- 
ations from rtfae grand canonical ensemble results occur 
at small Ao.Ea The important point, though, is that the 
zeros of the harmonics in the two schemes agree. 



IV. ABSENCE OF PERSISTENT MAGNETIC 
OSCILLATIONS 

The calculations in the previous section discussed the 
variation of different properties of the vortex lattice state 
within one particular period (n M £ (20,21)) of the nor- 
mal state magnetic oscillations. In order for the mag- 
netic oscillations to persist in the vortex lattice state, the 
same pattern of variation must occur over many periods 
of the normal state magnetic oscillations. To investigate 
whether or not this is the case we have studied the de- 
pendence of superconducting properties on field through 
a number of periods of the normal state oscillations. The 
small Ao behavior always involves the quantity Aq /nj/ 4 
and retains the normal state magnetic oscillations with 
reduced amplitude. The zero and subsequent sign re- 
versal of the fundamental harmonic with increasing Ao, 
however, does not occur in every normal state oscillation 
period. In Fig. 5, we plot the sine of the fundamen- 
tal harmonic versus for a value of AA equal to 4.75 
(where the weak maximum in Fig. 4 occurs). These re- 
sults show that no clear component of the magnetization 
with the normal state period survives in the crossover 
regime. The harmonic of the Fourier tranform of the 
magnetization in the finite interval from to n M + 1 in 
this regime varies in sign and magnitude with no pat- 
tern we have been able to discern, consistent with results 
presented in our earlier work.Ej Corresponding variations 
occur in the lowest band quasiparticle energies. In Fig. 

6 we plot the eigenvalue means as in Fig. 3 but for the 
case n M £ (24,25). One sees that the three curves con- 
verge together as in Fig. 3 but this time do not cross 
when the crossover regime is entered. That is, the cross- 
ing effect of Fig. 3 may or may not occur depending on 
Landau level index. To emphasize this, we plot in Fig. 

7 the fundamental harmonic averaged over two different 
six period intervals. We see that the phase shift effect 
of Fig. 4 has been completely washed out and the fun- 
damental harmonic is smoothly damped to zero as Fq 
increases beyond Tiuj c . 

V. CONCLUDING REMARKS 

Experimental evaluations of dHvA oscillation ampli- 
tudes are based on Fourier transforms over many periods 
of oscillations. The results in the proceeding section in- 
dicate that no measureable oscillation with the normal 
state period or with any other period we have been able 
to recognize occurs once the typical value of the pair- 
ing self-energy becomes comparable to the Landau level 
separation. Because we work with relatively small Lan- 
dau level indices compared with the typical experimen- 
tal situation, we are not able to completely eliminate 
the possibility that oscillations in this regime are peri- 
odic with-a different periodicity or with a periodicity in 
B ratherO than in B^ 1 , although we have looked for 
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such patterns without success and are reasonably confi- 
dent that they do not exist. It seems clear that in the 
3D case where many Landau levels contribute even for a 
fixed field, that magnetic oscillations will be even more 
strongly suppressed. Disorder broadening, which we have 
neglected, will damp the oscillatory signal beyond that 
calculated here. 

In conclusion, we have done a detailed analysis of the 
nature of the quasiparticle states in the field regime near 
the upper critical field of a 2D type-II superconductor. 
We find that for small A, all Landau levels and not just 
the Landau level at the Fermi energy are split. This 
property is associated with the absence of time-reversal 
symmetry in the presence of a magnetic field. The split- 
ting would be naively expected to lead to a sign change 
in the fundamental harmonic of the Fourier transform of 
the magnetization for a value of the pairing self-energy of 
order the cyclotron energy, analogous to the sign changes 
which can occur due to the spin-splitting of Landau lev- 
els. However our numerical calculations show that once 
the pairing self-energy is comparable to the normal state 
Landau level separation, although the spectrum of quasi- 
particle excitations and the magnetization have sizable 
variations on the magnetic field scale of the normal state 
dHvA oscillations, the variations are aperiodic. Accord- 
ingly, we find that dHvA oscillations are strongly sup- 
pressed once this regime is reached. 
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FIG. 1. Lowest eigenvalue (hu c units) vs AAo in an ap- 
proximation where all pairing matrix elements are the same 
constant [— l/2(n M 7r) _1,/4 Ao7kc>c] for n M = 20 (solid points), 
n M = 20.25 (pluses), and n M = 20.5 (open points). 



FIG. 2. Density of states versus energy (hu> c units) for 
AAo = 1 and n M = 20. 

FIG. 3. Mean eigenvalue (%uj c units) of the lowest band vs 
AAo for the flux lattice. Same notation as Fig. 1. 



G 



FIG. 4. Fundamental sine harmonic of M(B) (solid points) 
and of the magnetization of -1/2 the mean of the first quasi- 
particle band (open points) vs AAo (grand canonical) from 
the period n M £ (20,21). The pluses are results for M(B) 
generated in the canonical ensemble. 

FIG. 5. Fundamental sine harmonic of M(B) vs n M for 
AAo =4.75. Each point represents a calculation over a single 
period. 

FIG. 6. Mean eigenvalue (Juo c units) of the lowest band 
vs AAo for n M = 24 (solid points), n M = 24.25 (pluses), and 
n M = 24.5 (open points). 

FIG. 7. Fundamental sine harmonic of M(B) vs AAo for 
the periods n M £ (20, 26) (solid points) and n M € (21, 27) 
(open points). 
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